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ANALOG TO DIGITAL CONVERSION SHIFT ERROR CORRECTION 

DESCRIPTION 

The present application relates to the digital 
processing arts. The invention finds particular application 
in connection with a technique for correction of 
inconsistencies arising during pulse integration in analog to 
5 digital converters in nuclear medicine due to shifts in the 
timing of samples and will be described with particular 
reference thereto. It will be appreciated, however, that the 
invention finds application in a variety of applications 
where an analog pulse is digitally integrated. 

10 Nuclear medicine imaging involves the detection of 

radiation, such as gamma photons, from a subject under 
investigation. Gamma photon detectors for use in nuclear 
medicine imaging typically include a number of 
photomultiplier tubes coupled to a scintillation, crystal. 

15 When a gamma photon interacts with the scintillation crystal, 
a cascade of lower energy photons, typically in the UV band, 
is emitted, termed an "event." The nearest photomultiplier 
tubes, in turn, detect the cascade of lower energy photons 
and each produces an analogous electrical signal in the form 

20 of an analog event pulse. A measure of the total energy of 
the pulse is then used, in combination with that of other 
pulses, to recreate an image depicting radio isotope 
distribution in the subject. 

The total energy of the pulse is related to its 

25 integral, i.e., the area under a plot of energy vs. time. In 
conventional analog systems, the energy of the pulse is 
stored as an analog change level which is then measured 
directly. More recently, analog to digital converters (ADCs) 
have been used to determine the integral. Rather than 

30 measuring the entire output as an analog event, the amplitude 

■ 

of the electrical pulse signal is digitized at some 
continuous sampling frequency, and the digital samples 
scanned digitally to generate the integrated value. When an 
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event takes place (as indicated by a separate timing 
channel), the digital values are summed for a fixed number of 
samples . This summation forms a representation of the 
integrated electrical signal. 
5 One difficulty with this approach is commonly 

referred to as ADC shift error. This arises from the 
inability to match the start of each pulse accurately to the 
sampling intervals. The time distribution of gamma photons 
interacting with the detector is continuous. The ADC 

10 samples, however, at a fixed sampling frequency. As a 
result, the ADC does not sample the electrical signal in a 
manner that is consistent from pulse to pulse. For a given 
pulse, the ADC samples can be shifted temporally by ± H of 
the ADC sample period. This is the ADC shift error. 

15 If the pulse has a pure Gaussian shape, and the 

total integration period encompasses the entire event period, 
then the ADC shift error does not result in a variation of 
the integral. In practice, however, the pulses have a 
Poisson distribution (generally, a sharp rise followed by a 

20 ' slower decline) . While the pulses are often reshaped in the 
analog domain prior to conversion, they still do not. follow a . 

i 

strictly Gaussian distribution. Because of the non-ideal 
pulse shape, the ADC shift error results in a variation in 
the determined energy of the pulse . 

25 The distribution of the variation in determined 

pulse energy adds in quadrature with the underlying variation 
of the scintillator and photomultiplier tubes. This 
effectively widens the observed distribution of event 
energies. Since event energy is a primary means of 

30 discriminating good events from bad events, a widening of the 
observed distribution of good events leads to a widening of 
the discrimination window, thus increasing the likelihood of 
contamination by bad events. If the discrimination window is 
not widened, a loss in sensitivity results. In either case, 

35 image quality is reduced. 

In conventional gamma cameras, the contributions of 
the scintillation crystal and photomultiplier tubes to the 
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energy resolution of the camera are greater than that of the 
ACD shift error, so the ADC shift error does not have a 
significant effect on the overall resolution. However, as 
the energy resolution of these components improves, the ADC 
5 shift error is making a more significant contribution to the 
overall resolution. 

The amount of variation in the energy distribution 
(energy resolution) is typically identified by the full width 
at half maximum (FWHM) of the distribution. For a • normal 

10 distribution, FWHM relates to the standard deviation of the 
distribution and is generally calculated by multiplying the 
standard deviation by 2x (2xln (2) ) 1/2 (approximately 2.35). 
Roughly, it corresponds to the amount of additional error in 
the pulse, above that of a theoretical pulse with no 
•15 intrinsic variation (a pure spike), by virtue of the ACD- 
shift error. The extent to which the ADC shift error 
increases the FWHM is determined by the pulse shape and the 
sample frequency of the ADC. For example, for a given pulse 
shape, an integration period (the time between samples) of 20 

20; ns and 5 samples, a FWHM of about 25% may be obtained. This, 
drops to 13% when the integration period is reduced to 16 ns 
and to 5% when an integration period of 10 ns is used. 
Further reductions to about 1.5% are achieved by increasing 
the number of samples from 5 to 8, at the 10 ns integration 

25 period. 

While shortening the integration period (reduces 
the FWHM, some amount of loss will always be present. 
Additionally, there are practical limits to the frequency of 
sampling. 

30 The present invention provides a new and improved 

method and apparatus which overcome the above-referenced 
problems and others. 

35 In accordance with one aspect of the present 

invention, a method of reducing errors resulting from a 
temporal shift between an analog pulse and digital sampling 
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intervals is provided. The method includes digitally 
sampling the analog pulse at a plurality of spaced sampling 
intervals to generate a set of digital samples. An integral 
is determined from the samples in the set. A subset of the 
5 samples is selected. A correction factor is determined from 
the subset of samples. The correction factor is applied to 
the integral to generate a corrected integral value . 

In accordance with another aspect of the invention, 
a system for reducing temporal shift errors between an analog 

10 pulse and a digital sampling interval is provided. The 
system includes a means for digitally sampling the analog 
pulse at sampling intervals to generate a set of digital- 
samples, a means for determining an integral of the analog 
pulse from the set of digital samples , a means for selecting 

15 a subset of the digital values , a means for determining a 
correction factor from the subset of digital samples, and a 
means for applying the correction to the integral . 

In accordance with another aspect of the invention, 
a system for correcting a signal is provided. The system 

20 includes means for sampling the signal at intervals of time 
to generate a set of test samples. A correction table 
includes a plurality of codes, each of the codes 
corresponding to a relationship between samples in a subset • 
of predicted samples. Each of the subsets of predicted 

25 samples is selected from a different set of predicted 
samples. Each of the sets of predicted samples is shifted in 
time by a time interval which is less than an interval 
between samples. Each of the codes is associated with a 
corresponding correction factor. Means are provided for (a) 

30 assigning a. code to the set of test samples corresponding to 
a relationship between a subset of the test samples, (b) 
accessing the correction table to determine a correction 
factor corresponding to the code, and (c) applying the 
correction factor to an integration of the set of test 

35 samples. 

In accordance with another aspect of the invention, 
a system for reducing temporal shift errors between an analog 
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pulse and a digital sampling interval is provided. The 
system includes an analog to digital converter for sampling 
the pulse at intervals of time to generate a set of digital 
samples. A correction table assigns a correction factor to 
5 each of a plurality of codes, each of the codes corresponding 
to a relationship amongst samples in a subset of calibration 
samples. Each of the subsets of calibration samples is 
selected from a different set of calibration samples. Each of 
the sets of calibration samples are shifted in time relative 

10 to a sampling interval. A processor calculates an integral of- 
the pulse from the digital values and the correction factor. 

In accordance with another aspect of the invention, 
a method for correcting a signal is provided. The method 
includes sampling the signal at intervals of time to generate 

15 a set of test samples. A correction table which includes a 
plurality of codes is accessed. Each of the codes 
corresponds to a relationship between samples in a subset of 
predicted samples. Each of the subsets of predicted samples ■ 
is selected from a different set of predicted samples. Each 

20 of the sets of predicted samples is shifted in time, relative 
to other sets of predicted samples, by a time interval which . 
is less than the interval between samples. Each of the codes 
is associated with a corresponding correction factor. A code 
is assigned to the .set of test samples corresponding to a 

25 relationship between a subset of the test samples. The 
correction table is accessed to determine a correction factor 
corresponding to the code. An integration of the set of test 
samples is performed and the correction factor applied to the 
integration of the set of test samples. 

30 One advantage of at least one embodiment of the 

present invention is that the effect of ACD shift error on 
pulse integration is reduced. 

Another advantage of at least one embodiment of the 
present invention is that it enables a time stamp for a pulse 

35 to be created. 
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Another advantage of at least one embodiment of the 
present invention is that energy resolution of a gamma camera 
is improved. 

Still further advantages of the present invention 
5 will become apparent to those of ordinary skill in the art 
upon reading and understanding the following detailed 
description of the preferred embodiments. 



10 The invention may take form in various components 

and arrangements of components, and in various steps and 
arrangements of steps . The drawings are only for purposes of 
illustrating a preferred embodiment and are not to be 
construed as limiting the invention. 

15 FIGURE 1 is a schematic diagram of a gamma camera 

and energy processing system according to the present 
invention; 

FIGURE 2 is an enlarged sectional view of one of 
the detector heads of FIGURE 1; and, 
20 FIGURE 3 is .an exemplary plot of a . normalized 

scintillation pulse, before and after Gaussian shaping. 

With reference to FIGURES 1 and 2, a nuclear 
25 medicine imaging device A, such as a gamma camera, includes 
one or more detector heads 10 (three in the embodiment 
illustrated in FIGURE 1) . Each head detects radiation, such 
as gamma photons, from a subject in a region 12 under 
investigation, such as a patient. Typically, the radiation 
30 is the result of the decay of a radiopharmaceutical, which 
decay emits a radiation photon of characteristic energy. 
Each of the detector heads 10 includes a number of 
photomultiplier tubes 14 optically coupled to a scintillation 
crystal 16. A gamma photon interacts with the scintillation 
35 crystal generating a pulse of light or scintillation. The 
event is recognized by a timing channel of an analog to 
digital converter (ADC) 18. Each photomultiplier tube 14 



6- 



WO 2005/069040 



PCT/IB2005/050043 



which views the scintillation produces an analogous 
electrical signal in the form of an event pulse which is 
proportional to the amount of received light. The ADC 18 
includes sampling software which samples the pulse at a 
5 selected frequency. For example, as illustrated in FIGURE 3, 
energy samples Si, S2, S 3 / . • . S 5 of the pulse at sequential 
times ti, t2, t 3 , . . . tio are taken. The ADC does not 
detect the exact starting point t 0 of the pulse (i.e., at 
zero time) . Thus, there is an ADC shift error of ti-to which 
10 has a value of between 0 (in the rare event that the first 
sample coincides exactly with the start of the pulse) and 
(t2~ti) . Shaping of the pulse is optionally performed prior 
to sampling of the pulse to create a more Gaussian 
distribution . 

15 A processor 22 receives the sampled digital pulse 

energy values, Si, S 2 , S 3 , . . . S 5 from the ADC. The 
processor performs an integration of the samples to obtain an 
uncorrected value for the total energy of the pulse e.g., Si 
+ S2 + S3 + S 4 + S5. The processor uses a subset of these 

20 samples to create a code which correlates with a relationship 
between the samples in the subset. For example, a . plurality 
of selected consecutive samples e.g., Si, S2, and S 3 ■ is 
selected to form the subset. The samples in the subset are 
preferably selected to encompass a portion of ■ the energy 

25 distribution in which the energy is changing rapidly, 
preferably spanning the peak energy P e . The processor 22 
accesses a correction algorithm table 24 to obtain a 
correction factor related to the code. The correction factor 
is then applied to the integrated value of the energy (e.g., 

30 the two are multiplied) to generate a corrected energy value. 
The effect of ACD shift error on the resulting integrated 
energies is thereby eliminated or reduced. 

The corrected energy values of all tubes who saw 
the same event are summed and the sum is sent to an image 

35 processor and memory 26, which compares the total energy of 
each scintillation with the characteristic energy of the 
radiopharmaceutical to identify valid signals. The spatial 
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location of the head, each viewing photomultiplier tube, and 
the relative amplitude P e of the signals from viewing 
photomultiplier for each valid event are reconstructured into 
a three dimensional image. A video processor 28 withdraws 
5 selected portions of the data from the image memory 26 to 
generate corresponding human-readable displays on a video 
monitor 30. Typical displays include reprojections, selected 
slices or planes, and the like. 

The correction factor is a function of a predicted 

10 relationship between the subset of samples (Si, S2, S3) of the 
input pulse, the shift error, and its effect on the 
uncorrected energy value. This predicted relationship is 
based on an understanding of the expected shape of the pulse . • 
Knowing the general shape of the signal , these relationships 

15 can be determined. As the ADC shift error is varied, the 
relationship changes in a manner which is calculable and 
which can be encoded. Additionally, the amount of variation 
in the energy signal can be estimated, relative to the energy 
signal at an ADC shift error of zero. A ratio of the 

20 estimated energy signal at an ADC shift error of j to an. 
estimated energy . signal at an ADC shift error of 0 (zero) ■ 
forms a set of correction factors. These can be selected and 
applied, based on the encoded relationship between the subset 
of samples in an actual pulse. 

25 To generate the algorithm or look-up table 24, a 

set of correction factors is determined from a predicted 
pulse and the corresponding relationship between a selected 
set of samples is encoded. 

The correction factors are determined from a 

30 predicted energy distribution curve of a pulse. The 
predicted pulse is one which is anticipated to be similar in 
shape to an actual pulse to be integrated. The predicted 
pulse distribution curve may be based on a priori knowledge 
of the pulse shape. In one method, a gated integrator is 

35 employed with an analog system to create a plot of an actual 
energy pulse from the nuclear medicine device A or a similar 
device. The gated integrator allows a pulse to be analyzed 
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with is little or no sample error. Other methods of creating 
a predicted energy distribution curve are also contemplated, 
including sampling a pulse with an ADC at very short 
frequencies, or trial and error methods. In one trial and 
5 error method, a proposed plot is created, e.g., a Poisson or 
Gaussian distribution. A set of correction factors is 
determined from the curve. The correction factors are applied 
to a plurality of real pulses. It is then determined whether 
an improvement in energy resolution is achieved. This 
10 procedure may go through many iterations, each iteration 
evaluated to determine whether it improves the energy 
resolution. 

In one embodiment, the correction factors and 
associated values for each of a plurality of shift errors are 

15 determined for the predicted pulse, integration period, and 
total number of samples, as follows: 

1. A set of energy samples (e.g., Si, S 2 , S 3 , . . 
. Sio) of the predicted pulse at times (ti, t 2 , t3, . . . tio) 
are taken for a plurality of shift errors j from 0 to the 

20 integration period p=(t2-ti). These samples may be taken 
from the raw pulse energies or, as illustrated in FIGURE 3, 
after smoothing of the curve to achieve a more Gaussian 
distribution. The shift errors are preferably equally spaced 
in time over the range of possible shift errors. For 

25 example, if the integration period p is J nanoseconds (ns) , 
correction factors for shift errors of 1, 2, 3, . . .J ns 
are obtained. The more closely spaced in time the shift 
errors, the more accurate the correction will tend to be. 
However, for practical purposes, determining correction 

30 factors for- about 6-25 spaced shift errors is generally 
sufficient. For example, for sampling at 10 ns intervals, 
selecting shift errors of 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, and 
10 ns is sufficient to provide a substantial improvement in 
FWHM, i.e., the energy resolution of the nuclear medicine 

35 device. The integration period p is the same for both the 
determination of correction factors and for the sampling of 
an actual pulse. 
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2. The predicted curve is integrated analytically 
at each of the shift errors to obtain energy values E 0 to E p . 

3. A subset of N samples is selected from the 
samples for each integration. For example, sequential 

5 samples Si, S 2 , S 3 of samples Si, S 2 , S 3 . . . Sio are selected 
(N=3) . The selected samples in the subset are the same as 
those which will later be used for the correction of an 
actual pulse. The samples in the subset may be consecutive 
or non-consecutive. The samples in the subset are preferably 

10 in a region of rapidly changing energy, spanning or adjacent 
to the peak. While all samples in the set could, in theory, 
be selected (and thus, for the purposes of this application 
be considered a subset) , it is not necessary to use all 
samples, and could result in an unmanageable amount of data. 

15 for the processor to handle. In general, about two to four 
samples in the subset is generally sufficient. 

4. Each sample in the subset is normalized to the 
maximum sample in the integration. For example, the highest 
energy is given a value of 1 and the samples in the subset 

20 are assigned a value of the ratio of the energy of ■ the sample 
to the energy of the highest sample. 

5. The value of each normalized sample is encoded 
as an M-bit value, where M is preferably between 3 and - 6. 
The value of M will depend, in part, oh the data handling 

25 capacity of the processor. A value of 4 provides an adequate 
range of code values while being amenable to current 
processor capabilities. For example, the values are all 
multiplied by a common factor (15 in the case of a 4 -bit 
code), and expressed as an integer. The integer is then 

30 converted to binary code. 

6. The N M-bit code values from each subset are 
concatenated to form an NxM bit code Cj (3x4=12 bit code, for 
example) that represents the relationship between the sample 
subset for that integration at ADC shift error j . For 

35 example, if the normalized and factor multiplied code values 
of a subset are 0, 8, and 15, when converted to a binary, 
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twelve bit code, they are expressed as 0000 1000 1111 (i.e., 
equivalent to the number 143 in the decimal system) . 

7. Correction factors are assigned for each of 
the plurality of shift errors j between 0 and the integration 

5 period p=(t 2 -ti) . Each 'correction factor is determined from 
the ratio of the integrated energy of the pulse at the shift 
error to the integrated energy of the pulse at a shift error 
of 0, i.e., 

Ratio Rj= Ep/Ej 

10 The correction factor has a value of 1 at zero 

shift error. As the shift error varies, the value of the 
correction factor changes, but not generally in a linear 
fashion. 

8 . Using the codes from step 6 and the ratios 
15 from step 7, the correction table 24 is generated, as 

illustrated in Table 1. 
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TABLE 1 



Shift Error 


Code 


Correction 
factor 


0 


Co 


R 0 =l 


• 
• 


a 
a 
• 


a 
a 
a 


J 







9. Where codes C 0 to Cj do not fill out the 
5 entire NxM bit codespace, the correction table is completed 
by assigning correction factors based on the nearest defined 
code or interpolating correction factors from two or more 

■ * 

nearest defined codes. 

Once the correction table 24 has been generated, it 
10 can be fed into the processor, if it has been created 
separately- The correction table 24 can then be applied to 
actual pulses . 

Application of the correction table to an actual 
pulse is done on an event by event basis according to the- 
15 following steps. 

1. A set of energy samples Si, S2/ S 3 , . . . Si 0 
of the actual "test" pulse at times ti, t2, t3, . . . tio is 
taken. These samples may be taken from the raw pulse 
energies or, as illustrated in FIGURE 3, after smoothing of 

20 the curve to achieve a more Gaussian distribution. 

2. The maximum energy value of the set of energy 
samples from the test pulse is determined. 

3. From the set of energy samples, a subset of 
energy samples is selected. The same subset as was used to 

25 form the correction table is used, e.g., samples Si, S2, S3 of 
the test pulse at times t if t 2 , t 3 . 

4 . The subset sample values are normalized to the 
maximum value in the set of samples, as described for forming 
the correction table. 
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5, The subset is encoded as described in steps 5 
and 6 above for the correction table to generate a NxM-bit 
code Cj . 

6. Using the correction table , the generated code 
5 Cj is matched by the processor to the corresponding 

correction factor. Where the code is not listed in the 
correction table , the next lowest numerical code is selected 
and the corresponding correction factor for that code is read 
from the table* 

10 7 . The integrated result from the entire set of 

samples is multiplied by the correction factor to generate a 
corrected energy value. In the case where piling up occurs 
(two or more pulses overlap such that the trailing edge of 
one pulse overlaps the leading edge of the next) , some of the 

15 samples in the set may be corrected to account . for the 
effects of piling up prior to integration. 

The correction of the energy value in this way 
increases the resolution possible with the nuclear medicine 
device A by reducing one potential source of error. It may 

20 be used in conjunction with other error reducing methods , ; 
such as decreasing the integration period, for example, from 
about 20 ns, as currently used, to about 10 ns. 
Additionally, the technique allows scintillation crystals 
with higher energy resolution to be used without the ADC 

25 shift error overwhelming the resolution benefits gained. For 
■ example, in place of conventional sodium iodide as the 
scintillation crystal, which has an energy resolution of only 
about 10%, lanthanum bromide, with an energy resolution of 
about 3%, may be used as the scintillation crystal, improving 

30 the energy resolution of the nuclear medicine system. 

In addition to correcting the integration of the 
energy, the correction table is also optionally used to 
determine the start time of the pulse . Knowing the code 
value, the corresponding ACD shift can be read from the 

35 table. By subtraction of the ACD shift from the actual time 
at which the first sample was taken, the time at which the 
pulse commenced is determined. Creating a time stamp for the 
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pulse in this way finds application, for example, in positron 
emission tomography (PET) systems, to create a time to 
digital conversion . 

Without intending to limit the scope of the 
5 invention, the following theoretical specific embodiment 
demonstrates the application of the correction method. 

Consider that a pulse of the shape described in 
FIGURE 3, after Gaussian shaping, is sampled. An integration 
period of 10 ns with 5 samples is used. Without application 
10 of the correction table, this results in a FWHM of 5%. 

To create the correction table, the pulse is 
analytically integrated with ADC shift errors of 0-10 nS, as 
shown in Table 2 for each of the five samples. The integral, 
Ej, is the sum of the energies for the five samples, at the 
15 particular shift error. It can be seen that this value 
varies from about 2.22 to 2.35, depending on the shift error. 
The maximum energy value in the set is also identified. 

Table 2 

20 



ADC 
Shift 
error , j 


Integral 


Si 


s 2 


s 3 


S 4 


S 5 


0 


2.2163242 


0 


0.020316 


0.8657-68 


0.894791 


0.435449 


1 


2.2216111 


0 


0.047723 


0.930051 


0.851284 


0.382553 


2 


2.294135 


0 


0.094291 


0.977576 


0.80574 


0.351807 


3 


2.2442245 


0 


0.0162645 


1.008819 


0.75925 


0.31351 


4 


2.2669492 


0 


0.0251724 


1.024902 


0.712473 


0.277851 


5 


2.2945707 


0 


0.356621 


1.027339 


0.665669 


0.244942 


6 


2.3213471 


2.04E-05 


0.0469779 


1.01782 


0.618892 


0.214835 


7 


2.3412209 


0.000284 


0.0583147 


0.0998076 


0.572198 


0.187517 


8 


2.35017858 


0.001748 


0.689921 


0.969837 


0.525792 


0.162878 


9 


2.3477059 


0.006918 


0.785175 


0.934842 


0.0480043 


0.140728 


10 


2.337165 


0.020316 


0.865768 


0.894791 


0.0435449 


0.120841 



Subsets comprising the first three samples Si, S 2 , 
S3 are used from each ACD shift error to generate a code. 
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Specifically, the Si, S2, S3 values are normalized to the 
maximum value at each shift error. For example, the maximum 
value for the zero shift error set is S 4 = 0.894791. Thus, 
the normalized value of S 3 is 0.865768/0.894791 = 0.967564. 
5 Each of the normalized values is then multiplied by a common 
factor: 15 in the present example. Thus, for the example 
given, 15 x 0.967564 = 14.51. This is rounded down to the 
nearest integer = 14, which is equivalent to the 4 bit code 
1110. 

10 The codes for the three values are then 

concatenated to form a 12 bit code. In the example, the 
corresponding values for Si and S 2 in the subset are both 
zero, so the concatenated code is 000000001110, equivalent to 
a decimal code of 14. This is the code assigned to a shift 

15 error of 0. Corresponding codes for each of the other shift 
errors are then determined/ as shown in Table 3. The 
correction factors are determined from E 0 /Ej . In the first 
subset, for example, Ej = E 0 , so this value is 1. 

20 Table 3 



Shift 
Error , j 


s x 


s 2 


S 3 


Code 

¥ 


Correct:! on 
Factor 


0 


0 


0 


1110 (14) 


1110(14) 


1 


1 


0 


0 


1111 (15) 


1111 (15) 


0.99762 


2 


0 


0001 (1) 


1111 (15) 


00011111 (31) 


0.994129 


3 


0 


0010 (2) 


1111(15) 


00101111(47) 


0.987568 


4 


0 


0011 (3) 


1111(15) 


00111111(63) 


0.977668 


5 


0 


0101 (5) 


1111 (15) 


01011111(95) 


0.965899 


6 


0 


0110 (6) 


1111 (15) 


01101111 (111) 


0.954758 


7 


0 


1000 (8) 


1111(15) 


10001111 (143) 


0.946653 


8 


0 


1010 (10) 


1111 (15) 


10101111 (175) 


0.943046 


9 


0 


1100(12) 


1111 (15) 


11001111(207) 


0.944038 


10 


0 


1110 (14) 


1111 (15) 


11101111 (239) 


0.948296 



This provides the correction table 24, which in the 
illustrated embodiment, is represented by Table 4. 

-15- 



WO 2005/069040 



PCT/IB2005/050043 



Table 4 



Shift Error , j 


Code 


Correction 


0 


14 


1 


1 


15 


0.99762 


2 


31 


0.994129 


3 


47 


0.987568 


4 


63 


0.977668 


5 


95 


0.965899 


6 


111 


0.954758 


7 


143 


0.946653 


8 


175 


0.943046 


9 


207 


0.944038 


10 


239 


0.948296 



5 The correction table is then entered into the 

processor 22 of the nuclear medicine device A, to be used on 
actual test pulses which are sampled at the same integration 
period and number of samples. 

To apply the correction table, the set of samples 

10 from each event is examined to determine the maximum sample. 
Samples in subset Si, S 2 , S 3 are then normalized to the 
maximum sample and 4-bit codes are generated. The 4-bit codes 
are concatenated to a 12 bit code and this is used to index 
the correction table. The resulting correction factor is 

15 multiplied by the integration of the set of samples to form 
the corrected energy value. 

In the illustrated example, the value of Si is 
always zero. Thus, it does not affect the relationship 
between the three samples in the subset. Where this is the 

20 case, a different subset, such as S2, S3, S 4 could be used, 
such that each of the three values plays a part in 
determining the relationship. 

To test the process, the test pulse was 
analytically integrated in this manner a thousand times using 

25 ADC shift errors that are uniformly distributed between 0 and 
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the integration period (10 ns in this case) . A partial list 
of those integrations is shown in Table 5. 

The resulting FWHM for the corrected data is 0.35%. 
This is a substantial improvement over the 5% FWHM without 
5 correction. It will be noticed that the method works even 
for negative ADC shift errors, i.e., those where the first 
sample was taken prior to the actual start time of the pulse. 

The invention has been described in terms of 
correction of energy values in an energy pulse from a gamma 

10 photon detector used in nuclear medicine imaging. This 
includes Single Photon Emission Planar Imaging, Single Photon 
Emission Computed Tomography (SPECT) , and Positron Emission 
Tomography (PET) systems. It will be appreciated that energy 
pulses from other sources may be similarly corrected for 

15 errors due to difficulties in beginning the sampling a 
uniform time interval after the start of the pulse. The 
application may also extend generally to any situation where 
a continuously randomly occurring signal of limited duration 
and known shape is sampled by a fixed frequency ADC and where 

20 the integral of the sampled signal is desired. The invention 
finds wider application in treating data obtained from other 
sources where a variable is plotted against a determinate 
such as electrical current vs. time, and where the variable 
follows a generally predictable variation. vs. the, 

25 determinate. Specific examples include thermal energy vs. 
time, concentration changes vs. time, energy vs. wavelength, 
and the like. 

The invention has been described with reference to 
the preferred embodiment. Modifications and alterations will 
30 occur to others upon a reading and understanding of the 
preceding detailed description. It is intended that the 
invention be construed as including all such modifications 
and alterations insofar as' they come within the scope of the 
appended claims or the equivalents thereof. 

35 
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Table 5 



ADC 
bnir t 

error , 

• 

3 


Si 


s 2 


s 3 


s 4 


s 5 


Integration 
(Uncorrected) 


Code 


Correction 
factor 


Integration 
(Corrected) 


6 


0.000 


0.048 


0.930 


0.851 


0.393 


2.222 


27 


0.998 


2.216 


4 


0.007 


0.785 


0.935 


0.480 


0.352 


2.348 


244 


0.948 


2.226 


-3 


0.000 


0.094 


0.978 


0.806 


0.352 


2.229 


38 


0.994 


2.216 


-1 


0.000 


0.866 


1.025 


0.712 


0.278 


2.267 


73 


0.978 


2.216 


-1 


0.000 


0.866 


1.025 


0.712 


0.278 


2.267 


73 


0.978 


2.216 


-3 


0.000 


0.094 


0.978 


0.806 


0.352 


2.229 


38 


0.994 


2.216 


5 


0.020 


0.866 


0.895 


0.435 


0.121 


2.337 


330 


0.948 


2.216 


0 


0.000 


0.357 


1.027 


0.666 


0.245 


2.295 


98 


0.966 


2.216 


-3 


0.000 


0.094 


0.978 


0.806 


0.352 


2.229 


38 


0.994 


2. 216 


-3 


0.000 


0.094 


0.978 


0.806 


0.352 


2.229 


38 


0.994. 


2.216 


-4 


0.000 


0.048 


0.930 


0.851 


0.393 


2.222 


27 


0.998 


2.216 


-2 


0.000 


0.163 


1.009 


0.759 


0.314 


2.244 


53 


0.998 


2.216 


2 


0.000 


0.583 


0.998 


0.572 


0.188 


2.341 


155 


0.947 


2.216 


-3 


0.000 


0.094 


0.978 


0.806 


0.352 


2.229 


38 


0.994 


2.216 


3 


0.002 


0.690 


0.970 


0.526 


0.163 


2.350 


188 


0.943 


■ 

2.216 


4 


0.007 


0.785 


0.935 


0.480 


0.141 


2.348 


244 


0.948 


2.226 
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